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I propose a new paired-wavefunction with a parameter 
that continuously interpolates Irom the ID Jastrow-product 
to the Majumdar- Ghosh dimer-wavefunction appropriate for 
the frustrated Heisenberg S = 1/2 antiferromagnet. This spin 
paired-state constructed in Sz basis is an alternative to the 
well-known resonating-valence-bond basis state for describing 
the S = ground-state with no apparent long-range spin 
order. Some numerical evidences are presented. 



The simplest non-trivial single-particle wavefunction in 
the first quantized form is the plane wave, exp(iA:-r), cor- 
responding to a free particle in space. For the many-body 
system the simplest ground-state wavefunction can be 
constructed by either anti-symmetrizing or symmetrizing 
the single-particle plane waves depending on the statis- 
tics of the particles being considered. 

The next non-trivial step is to in-cooperate the ef- 
fects of interaction by taking (anti-)symmetric pairwise 
wavefunction and further (anti-)symmctrizing over all 
the other remaining variables in the many-body wave- 
function. The simplest such wavefunctions are the 
Jastrow-product [Q and the less well-known BCS paired- 
wavefunction projected onto the Hilbert space of defi- 
nite number of particles |g] . Recently, these two- types of 
wavefunctions have been applied to the fractional quan- 
tum Hall states H . In this letter I propose a new paired- 
wavefunction with a parameter that interpolates from 
the Jastrow-product which is an exact ground-state of 
the Haldane-Shastry model (HSM) [Q to the Majumdar- 
Ghosh wavefunction for the frustrated SU{2) Heisenberg 
spin chain ||^. 

One of the main results of this letter is the following 
wavefunction for the spin chain with Na sites 
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where M — Na/2 is the total number of down or up 
spins, the sum is over all the spin configuration a in Sz 
basis, and rrii (rii) the integer site index of ith down(up) 
spin. The second factor in the sum corresponds to the 
Marshall sign. The orbital function (p{z, a) is given by 



cj,{{z,a)=Bct[fiizl,zi,^)], 
fiix,y) = ix + yy/''\x-y\\ 
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where Zm = exp(i2TTm/Na) and I some integer. Here, 
Det denotes determinant of the M x M matrix [f{xi, yj)] 



while f{xi,yj) is a wavefunction for a pair of up- and 
down-spins. One can consider A > as a measure of 
apparent repulsion between the spin pair. 

I summarize some noteworthy features of the wave- 
function as follow: 

• The wavefunction is a SU(2) spin singlet and an 
eigenstate of the lattice translation operator with 
the wavenumber Q = 2-k{M — l)/Na- Hence, I = 
(/ — M) corresponds to Q = tt (Q = 0). 

• When A = Af - 1 and Z = for odd M (A = M-2 
and I = M for even M) the wavefunction exactly 
reduces to the following fully Gutzwiller-projected 
free-fermion wavefunction Ipll 
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which is the exact ground-state of the isotropic 
SU{2) Haldane-Shastry spin chain. The HSM is 
known to belong to the same universality class as 
the Heisenberg spin chain with the nearest neighbor 
exchange 



• When A = 1 it becomes the exact Majumdar-Ghosh 
wavefunction with definite crystal momentum Q ~ 
vr (for ; = 0) or Q = (for Z = M), 

l*> = l^>T|S), (5) 

where |^) = (Ti - it) • • • (Ti - it) and \B) = 
T\A). [T is the lattice translation operator.) 

• The wavefunction is in the Dyson's form of the 
BCS- wavefunction projected onto the Hilbert space 
of definite particle number. One major difference 
here is that the power A is positive, reflecting ap- 
parent repulsive interaction between the spin pair. 

One of the most important features of the proposed 
wavefunction mentioned above is the spin singlet prop- 
erty which is also a general requirement for the ground- 
state of antiferromagnetically correlated spin systems [0 . 
The wavefunction in Eq. (Q) is by construction 5*2 = 
state (i.e. M = Na/2); thus, if 5+|*) = it is also a 
S — Q state. This feature can be shown by proving that 
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where P'^ is the spin exchange operator. Using the an- 
tisymmetric property and minor expansion of the deter- 
minant it is straightforward to show that the singlet re- 
quirement reduces to the foUowing condition for some 
fixed focation zj, 

M 

X sgnO' - n)sgn{iy - i) f {z] , zj;^) f {zl , zj)Mii; , (7) 

where -A^I-^lv, Y is a (M — p) x {M — p) matrix ob- 
tained from the original M x M matrix by eliminating 
the columns ii,i2, ■ ■ ■ ,ip and the rows ji , J2, ■ • • , jp- This 
condition holds true for any even function / since the 
second term in Eq. (0) vanishes identically when summed 
over the dummy indices j and /i and the first term be- 
comes just a simple minor expansion of the original deter- 
minant. Thus, the wavefunction in Eq. (|l|) with / given 
by Eq. (g) is a particular case of SU{2) singlet. It is in- 
teresting to note that the singlet constraint is naturally 
resolved with the BCS-like form of wavefunction in the 
Sz basis, an alternative to the resonating-valence-bond 
basis. 

The pair function / with A > in the wavefunction 
suggests that the interaction between up- and down-spins 
is of repulsive nature. But, it is somewhat misleading 
since for the lattice system there could be an arbitrary 
constant factor ni<,(^i ~ ^j)™ ^P ^ ^ correction to the 
phase factor. Thus, even for the apparently repulsive 
case of A = Af — 1 (or A = M — 2), the wavefunction 
can be divided by the constant factor with m — 1 and be 
rewritten as 
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where the effective attraction between up- and down- 
spins (i.e. the short-range antiferromagnetic correlation) 
is made explicit. It is remarkable that for the case of 
A = 1 there also exists an equivalent alternative form 
given by 
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where A/" denotes proper normalization. In this repre- 
sentation the attraction for a spin pair appears to be 
infinite. While the first state corresponds to the ground- 
state of the gapless spin liquid where the spinous, the 
elementary excitations, are unbound and can have arbi- 
trarily low energy (i.e. the effective attraction in this 
case is not enough to subdue the pair-breaking quantum 
fluctuation and to induce a gap) , the second corresponds 
to the fully paired or dimerized state with gap in the 
energy spectrum. I can, thus, infer from the two known 
cases that as A decreases the effective attraction increases 



and that the Majumdar-Ghosh limit (A — 1) can be 
viewed as the strongest-possible-pair condensation. In 
fact the first case can be viewed as the critical point of the 
frustrated Heisenberg antiferromagnet given in Eq. dTfl) 
1^,0, and any smaller A corresponds to the massive 
dimerized phase. This transition can be thought of as 
due to the singlet pair condensation and that the size of 
the singlet pair fluctuation becomes smaller and smaller 
as the frustration is increased towards the Majumdar- 
Ghosh limit. 

In order to show more clearly the nature of pairing in 
the frustrated Heisenberg antiferromagnet I start with 
the following spin Hamiltonian 
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where S — 1/2 and 7 > is the frustration parameter. 
The Hamiltonian with 7 = corresponds to the Bethe- 
ansatz solvable model H]. At 7 = 1/2 the exact dou- 
bly degenerate ground-states are the Majumdar-Ghosh 
wavefunction [g . When 7 is greater than a critical value 
(« 0.241) the otherwise marginally irrelevant umklapp 
processes become relevant perturbation to the gapless 
spin-liquid state and the spin-chain spontaneously dimer- 
izes; and, the critical point is identified with the /?^ —^ 8tt 
limit of the sine-Gordon model Q . Employing the stan- 
dard method |l^,P one can reduce the Hamiltonian to 
the exactly solvable Luttinger model plus the following 
umklapp terms expressed with the Jordan- Wigner lattice 
fermions, 



nS=l,2 



n5=±l 

Hernfitian conjugate. 
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where oi = 1 and 02 — —7. The rapidly oscillating 
modes at the Fermi points have already been factored out 
from the fermion operators. When the z-axis anisotropy 
is allowed this umklapp term is also responsible for the 
Neel ordered state [ [12[ . 

Recently, Wilczek and Nayak interpreted the umklapp 
term in the two-dimensional Mott-system as pairing in- 
teraction of particle and hole and when the order pa- 
rameter (ipLipjj) is a pure imaginary number a BCS-like 
state rather than the charge-density-wave state is real- 
ized [Ol. In the spin language their BCS-like state cor- 
responds to the spontaneously dimerised state M, the 
charge-density- wave state to the Neel ordered state [O , 
and their particle-hole pairing to the spin singlet pair- 
ing. Furthermore, their BCS-like state through repulsive 
channel seems consistent with our proposed wavefunc- 
tion with A > which exhibits the apparent repulsive 
interaction. 



The original estimate of the critical frustration 7c = 
1/6 yi which can be obtained from Eq. (^l|) by Taylor 
expansion of the operators like 'ijj{x + S) is modified to 
7c — 0.277 if the full quantum operator-product expan- 
sion is used. (Numerical estimate of 7c is approximately 
0.241 [H.) For estimating 7c it is more convenient to 
change the dynamical variables to the following canoni- 
cally conjugate phase fields 



eix)=27T / dx'ipR{x')+pLix')), 



x)^TT dx'{pii{x') - pl{x')), 
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where pr.l are right- and left-density operators and 
[(j){x),9{x')] = i7rsgn(a:: — x'). The original spin SU{2) 
operators are represented in these fields as S'^{x) ex 
exp(±i(/)(x)) and Sz{x) ex exp(±z0(a;)). Furthermore, the 
lattice Jordan- Wigner fermions are decomposed into the 
right- and left-moving Mandelstam fermions represented 
in terms of the phase fields as follow 






ia4)(x) .. ij36(x) . 
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where :: means normal ordering and a/? = 1/2 for 
the fermionic statistics. In order to obtain the correct 
anticommutation relations between the right- and left- 
moving operators the zero-modes of the phase fields are 
separated as 6'(a;) = 9o{x) + 9{x) and ^(a;) = (f>o{x) + (j>{x) 
where 6*0 (x) = 6*0 + 2ttN/L, (j)o{x) = 0o + ^ttJ/L and 
[N, (j)o] = [J, 9o] = i Q- The operators 6*0, ^o, N, and J 
are associated with the bulk modes of the ID fluid and 
commute with the q ^ bosonic modes. (N here is not 
the number of particles but rather the fluctuations from 
the mean value A'o.) Some care in ordering the operators 
is needed when contracting the right- and left-moving 
operators and they are given as follow 



lim 'ip]^{x)ipLiy) = ysgn(x - y) (m^ + A^) 
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where r = 2/3^ — (5^/2, J3 = /3exp((^), a = Q;exp(— (/?) 
and ao is the lattice constant. A similar form appears in 
Ref. p6| but is missing the important factor sgn(a; — y) 
and the corrections due to the short-distance cut-off. The 
dimcnsionless parameter A is introduced to control the 
short-distance behavior for normal ordering two exponen- 
tial operators and, in this linear approximation, is fixed 
by the filling fraction and is given as A = 2/(e'-^7r) sa 0.357 
(C = the Euler's constant). If |77i| 3> A then the control 
parameter can be dropped from Eq. (HO) , but if the dif- 
ference is of the order of the lattice constant it cannot be 
ignored. 



In order to reproduce the correct anticommutation re- 
lations at the non-interacting value of the renormaliza- 
tion constant exp(— 2(p) = 1,1 have to set j3 = 1/2, a = 1 
and, thus, r = 0. More generally, I obtain the following 
anticommutation relations for the left and right Mandel- 
stam modes 

{ijRix),i^Uy)} = IM^). i^iiy)} = zs{x ~ y), (17) 

where Z = (r(i/ - l/2)/^T{i'))\^-^'' with i^ = (a/2 + 
/3)^. Thus, for the free fermions (i.e. exp{±ip) = 1) Z = 1 
or i^ = 1. And, this requirement along with the fermion 
constraint a/3 = 1/2 fix the values of a and /3. In the 
gapless phase of the spin chain, however, exp(— 2iy9) = 2 
]ll|,|. Thus, /32 = 1/8, a2 = 2 and T = -3/4. The umk- 
lapp terms Hu given in Eq. (|T^) can now be contracted 
to the following form |17t 



/1 + A2 
4 + A2 
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cos(26l(a;)) :, (18) 



where some overall constants have been dropped. The 
coefficient vanishes when 7 « 0.277 and changes sign 
if 7 further increases. The difference between this value 
and the numerical estimate is presumably due to the high 
energy modes with non-linear dispersion relation near the 
Brillouin zone boundary. 

For half-odd integer spin chains the Lieb-Schultz- 
Mattis (LSM) theorem dictates that the gap open for 
both of the two low-energy sectors at the crystal momen- 
tum and TT 11^ . And, if the ground-state is doubly de- 
generate one can construct two linear combinations that 
are not eigenstates of the lattice translation operator. 
Therefore, the dimer gap induced by the singlet spin- 
pairing should also accompany the spontaneously broken 
translational symmetry. In the presence of holes or in 
two-dimensional systems, however, the spin paired-state 
could exist without the broken lattice translational sym- 
metry. 

I now numerically demonstrate that the proposed 
wavefunction given in Eq. (Q) is a good trial wavefunc- 
tion. I do this by taking the overlap of the wavefunc- 
tion with the exact numerically found ground-state of 
Na = 16 spin chain given in Eq. dlG) with the frustration 
7 ranging from to 1/2 as shown in Table |. Listed in 
the table are the maximum overlap squared for the sec- 
tor Q = and Q = it for four different values of A and 
the corresponding optimum values of 7. (For M even 
(odd) and 7 < 7c the lowest-energy state in the sector 
Q = TT (Q = 0) is a 5 = 1 state and the 5' = state 
is the next lowest for any finite-size system. In all other 
cases of interest the ground-states for both sectors are 
spin singlets.) As shown in Table | the overlap squared 
is very close to one although it somewhat deteriorates in 
the region 7 w 0.4. 

In conclusion, I propose a new paired-wavefunction for 
the frustrated Heisenberg antiferromagnetic spin chain. 



The wavefunction continuously interpolates from the 
gapless spin liquid phase to the fully dimerized phase 
of Majumdar-Ghosh. From the bosonization method I 
estimate the critical value of the frustration and further 
argue that the gap is induced by the singlet pairing in 
accord with the proposed paired-wavefunction. I also 
give numerical evidences that show the validity of the 
paired-wavefunction as ground states of the frustrated 
spin chain. Details will be published elsewhere p7[ |. 

I further speculate that an analog of the trial wavefunc- 
tion given in Eq. (|l|) adopted to the 2D lattice (i.e. take 
Zj's as general complex numbers and use elliptic func- 
tions) is a possible ground-state for the 2D frustrated 
spin system. This S = construction in Sz basis is to 
be contrasted/compared with the RVB states used, for 
example, by Liang et al ||l9[. It is also tempting to con- 
jecture that the spin-spin singlet pairing correlation sur- 
vives the doping of the spin system with small density of 
holes and that the "built-in" superconductivity naturally 
arises in some reasonable parameter space. Numerical 
works are in progress. 
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TABLE L Square of the overlap of the trial wavefunctions 
given in the text with the exact 5* = ground-states of the 
spin-chain with Na ~ 16 sites at the wave-number Q = 
(or I = M) and Q = tt (or Z = 0). The frustration 7 = 1/2 
corresponds to the Majumdar-Ghosh point and 7 = to the 
Bethe-ansatz solvable Heisenberg spin chain. The overlaps 
for Q — n are slightly better than those for Q = in the 
intermediate region 0.25 < 7 < 0.5. 
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